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C^The effect of small forced symmetry breaking on the dynamics near a structurally stable heteroclinic 
* '"^ cvcle connecting two equilibria and a periodic orbit is investigated. This type of system is known 
£^to exhibit complicated, possibly chaotic dynamics including irregular switching of sign of various 

1 p'hase space variables, but details of the mechanisms underlying the complicated dynamics have 
jiot previously been investigated. We identify global bifurcations that induce the onset of chaotic 

^dynamics and switching near a heteroclinic cycle of this type, and by construction and analysis 
of approximate return maps, locate the global bifurcations in parameter space. We find there is a 



-^ threshold in the size of certain symmetry-breaking terms below which there can be no persistent 
f~y witching. Our results are illustrated by a numerical example. 

oo : 

o . 

\o . 

Introduction 

• *~{r f, is well-established that the presence of symmetries in dynamical systems 
^can result in the existence of heteroclinic cycles that are structurally stable 
.^with respect to symmetric perturbations [1,2]. By heteroclinic cycle we mean 
^a collection of two or more flow invariant sets {£i, . . . , of some system of 
C^ordinary differential equations together with a set of heteroclinic connections 

_ _ i 

{71 (t), . . . , 7n(t)}, where 7j(t) — > £j as t — > —00 and 7j(i) — ► as t — > +00, 
and where £ n +i = £i- I n many studies, all the £j are equilibria, but in this 
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paper we explicitly consider the case that one of the £j is a periodic orbit. The 
connections -fi may be isolated, or there may be a continuum of connections 
from £j to for one or more i. 

There is a large literature on structurally stable heteroclinic cycles (SSHC), 
including work establishing conditions for the existence and asymptotic stabil- 
ity of heteroclinic cycles [3-5], examination of the dynamics near heteroclinic 
cycles and networks of heteroclinic cycles [6-9] , and unfolding of bifurcations 
of heteroclinic cycles [10-12]. SSHC arise naturally in mathematical models of 
physical systems with symmetry or near-symmetry [13-16]. In these models, 
the physical system is idealised as having perfect symmetry, leading to the 
existence of invariant subspaces in the model and thus to the robustness of 
heteroclinic cycles with respect to symmetric perturbations. It is natural to 
ask how much of the dynamics observed in symmetric models persists under 
non-symmetric perturbations. Some effects of small symmetry-breaking have 
been documented [17-20], and aspects of the related question of how much 
of the dynamics persists under the inclusion of small noise have also been 
considered [21,22], but details are likely to vary greatly between different ex- 
amples. A few cases of experimental observation of near-heteroclinic cycles 
have been reported, most recently in [23], but see also the references therein. 
In these cases, experimental noise and small symmetry-breaking effects pre- 
vent exact heteroclinic cycles from occurring, but there is clear evidence for 
near-heteroclinic structures in certain regimes. 

Our interest in the particular set-up explored in this paper is motivated 
by [20], which makes the observation that the addition of small symmetry- 
breaking terms to a system containing a heteroclinic cycle connecting two 
equilibria and a periodic orbit (as well as symmetric copies of the cycle) results 
in seemingly chaotic dynamics, with orbits passing near the various equilibria 
in the system repeatedly but in an irregular pattern, as illustrated in figure 1. 
A main point of [20] was to show that repeated switching of orbits in this 
manner could arise in a simple four-dimensional, nearly symmetric model, 
but the specific mechanisms underlying the complicated dynamics were not 
explored in detail. 

In this paper, we examine a generalisation of the situation from [20], focusing 
on the structure and origin of chaotic dynamics in the system and on how 
switching dynamics is induced. Here and elsewhere in the paper, switching 
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Figure 1. Irregular switching in the time series of a dynamo model studied in [20, figure 1]. 
Panels (a) and (b) show the evolution of different coordinates of the same trajectory, and panel (c) 
shows the same coordinate as in (b) over a longer time interval. 

refers to the itinerary that an orbit follows under the dynamics. Specifically, 
in the fully symmetric version of our system there is a heteroclinic network 
consisting of four symmetric copies of the basic heteroclinic cycle. Invariance 
of various subspaces ensures that an orbit may make repeated passes near only 
one cycle. Once the symmetries are broken, however, an orbit may switch, i.e., 
make traversals near more than one of the original cycles (although, of course, 
the cycles themselves may not persist when the symmetry is broken). 

A main result of this paper is that in the case of small symmetry breaking, 
switching in one variable occurs when a complicated attractor arising from the 
presence of transverse homoclinic orbits of a periodic orbit crosses the stable 
manifold of one of the equilibria in the system. The existence of the transverse 
homoclinic orbits depends on a broken rotation symmetry, while the proxim- 
ity of the attractor to the stable manifold of the equilibrium is caused by a 
broken reflection symmetry. Switching in a second variable results from the 
interaction between broken reflection symmetry and complicated dynamics as- 
sociated with a heteroclinic bifurcation between the equilibria. Thus, switching 
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results from the right combination of a global bifurcation and small symmetry 
breaking. 

A second significant result of this paper is the observation that there is a 
threshold for the size of symmetry breaking below which persistent switching 
cannot occur. More precisely, the existence of the heteroclinic cycle requires 
three separate symmetries to allow structurally stable connections within three 
invariant subspaces. We control the degree to which the three symmetries are 
broken by three small parameters, e±, €2 and €3; ei controls the degree to which 
the periodic orbit in the cycle deviates from a perfect circle, while €2 and €3 
break reflection symmetries. For fixed small €2 and £3, we find that there is 
a threshold in t\ for persistent switching to occur. For sufficiently small ei, 
there may be a single switch from one part of phase space to another, but it 
is only for t\ beyond the threshold value that an orbit can repeatedly visit 
different parts of the phase space. We find that it is possible to get sustained 
switching in one or other or both of the variables associated with the reflection 
symmetries, and that the threshold values of ei are different for switching in 
the two variables. The threshold does not go to zero as €2 an d €3 go to zero. 

Sustained switching of orbits near heteroclinic cycles and networks has been 
observed in a number of other settings. Clune and Knobloch [24] describe 
an example in which there are two symmetrically related copies of a non- 
asymptotically stable heteroclinic cycle, with nearby orbits making repeated 
passes near each cycle; no mechanism for the switching is suggested in this pa- 
per. Aguiar et al. [8] find switching near a hybrid heteroclinic network formed 
from transverse heteroclinic connections between equilibria and connections 
that are robust because of symmetry; switching seems to result from the fold- 
ing and stretching caused by passage near the transversal heteroclinic connec- 
tions and by mixing near an equilibrium solution with complex eigenvalues. 
Kirk et al. [25] have an example of switching near a heteroclinic network that 
has no transversal connections; the switching is caused entirely by passage 
near an equilibrium with complex eigenvalues. Postlethwaite and Dawes [9] 
describe a variant of switching near a heteroclinic network in which each cycle 
in the network is unstable along a direction transverse to the cycle; orbits visit 
cycles in the network in a fixed order (being pushed away from each cycle in 
the transverse direction, which also happens to be the contracting direction 
for the next cycle) but the number of traversals of each cycle before switch- 
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ing to the next cycle can be constant or irregular. Ashwin et al. [26] describe 
switching associated with a stuck-on heteroclinic cycle between two invariant 
subspaces; here the switching is caused by a nonlinear mechanism that chooses 
between the different possibilities in a manner that is well modelled by a ran- 
dom process. Switching can also be induced by adding noise to a structurally 
stable heteroclinic network [22]; noise sensitive switching has been observed 
by [27,28]. None of these examples explicitly considers symmetry breaking as 
a mechanism for switching. 

We adopt a standard approach to analysis of the system of interest, i.e., 
we set up a simple symmetric model in which there exists a heteroclinic cycle 
connecting two equilibria and a periodic orbit (Section 2), construct a return 
map that approximates the dynamics near such a cycle, and then add generic 
symmetry breaking terms to the return map (Section 3, with details in the Ap- 
pendix). Analysis of the return map is fruitful in cases where partial symmetry 
is retained, and allows us to prove the existence and asymptotic stability of 
periodic orbits, quasiperiodic solutions or heteroclinic cycles in various cases 
(Sections 4.1-4.3). In the completely asymmetric case, the return map is in- 
tractable, but we are able to make predictions about the dynamics by assuming 
there is a generic unfolding of the partially symmetric cases (Section 4.4). The 
example discussed in Section 5 confirms and illustrates the analysis. Some 
conclusions are presented in Section 6. 

A complicating factor in the analysis presented in this paper is that the 
unstable manifolds of one pair of equilibria and of the periodic orbit are two- 
dimensional, and there are continua of heteroclinic connections along some 
parts of the cycle in the fully symmetric case. Linearising about a single hete- 
roclinic connection is not appropriate, and the usual method of analysis needs 
to be adapted to keep track of orbits in a neighbourhood of all the connections. 
Our approach is similar to that taken in [7,25,29]. We note that our analysis 
need not consider the issue of which connection from a continuum is selected 
by the dynamics (as investigated in, for instance, [7,26,30,31]) since in our 
case breaking of the symmetries forces a discrete set of transversal connections 
to be selected from each continuum. Note also that some results about the 
dynamics near a heteroclinic cycle connecting an equilibrium and a periodic 
orbit in a generic (i.e., non-symmetric) setting are described in [32,33], but 
the phenomena described in those papers will not be seen for small symmetry 
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breaking in our setting, and is not the focus of our interest here. 
2 Description of the problem 

We consider a system of ordinary differential equations x = f(x) where f : 
1R 4 — > IR 4 , and x = (x±, y±, x 2 , £3) G K 4 . It is sometimes convenient to use 
polar coordinates (ri, B\) such that z\ = xi + iyi = rie ldl . Initially, we assume 
the system is equivariant with respect to the action of a rotation and two 
reflections: «j(f(x)) = f(/Cj(x)), i = 1,2,3, where 

k 1 : (z 1 ,x 2 ,x 3 )^(z 1 e l4 ',x 2 ,x 3 ), 
k 2 : (z 1 ,x 2 ,x 3 )^(z 1 ,-x 2 ,x 3 ), 
k 3 : (z 1 ,x 2 ,x 3 )^(z 1 ,x 2 ,-x 3 ), 

with < (j) < 27T. These symmetries generate the group S 1 x Z 2 x Z 2 , and 
their presence ensures the existence of some dynamically invariant subspaces. 
We make the following assumptions about the dynamics in the subspaces, as 
illustrated in figure 2: 

• There exists a hyperbolic periodic orbit P in the invariant plane x 2 = x 3 = 0. 
Within this plane, the periodic orbit is a sink. 

• There exist hyperbolic, symmetry-related pairs of equilibria ±£^2 and ±E 3 
on the invariant lines z\ = 0, x 3 = and z\ = 0, x 2 = respectively. Within 
these lines, the equilibria are sinks. 

• Within the invariant subspace X3 = 0, P is a saddle and ±E 2 are sinks, and 
there are two-dimensional manifolds of heteroclinic connections from P to 
±E 2 (figure 2a). 

• Within the invariant subspace z\ = 0, ±E 2 are saddles and ±E 3 are sinks, 
and there are one-dimensional heteroclinic connections from +E 2 to ±E 3 , 
and from —E 2 to ±E 3 (figure 2b). 

• Within the invariant subspace x 2 = 0, ±E 3 are saddles and P is a sink, and 
there are two-dimensional manifolds of heteroclinic connections from ±£3 
to P (figure 2c). 

In the presence of the rotation symmetry k±, the coordinate 9\ decouples 
from the other coordinates, leaving an equivalent three-dimensional system 
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(a) (b) (c) 




Figure 2. The heteroclinic cycle for the fully symmetric system, (a) One of the connections in the 
x:i = subspace, from the periodic orbit P to the equilibrium point +E2; (b) the single connection 
in the x\ = yi = subspace, between the equilibria +E2 and +-B3; (c) one of the connections in 

the X2 = subspace, from +E3 to P. 



containing a SSHC connecting three saddle-type equilibria. This cycle may be 
asymptotically stable, depending on the eigenvalues at the three equilibria [4]. 
The behaviour of trajectories near such a heteroclinic cycle is well understood, 
with a typical orbit passing near each of the equilibria in a cyclic manner, 
spending ever increasing periods of time near each equilibrium. The dynam- 
ics in the fully symmetric, four-dimensional problem therefore has analogous 
behaviour: trajectories cycle between two equilibria and a periodic orbit, with 
the time spent near each equilibrium or the periodic orbit increasing with each 
subsequent traversal of the cycle [20]. Moving to four dimensions does more 
than replace one pair of equilibria by a periodic orbit: it also introduces dynam- 
ical features that will be important once symmetry is broken. In particular, as 
can be seen in figure 2, ±£2 and ±£3 are saddle-foci in the four-dimensional 
problem, and P and ±£3 have two-dimensional unstable manifolds. 

A detailed analysis of the effect on the dynamics of small symmetry break- 
ing is performed in the following sections; here we describe some geometric 
effects. Since ±£2, ±£3 and P are assumed to be hyperbolic in the fully sym- 
metric case, they persist and are hyperbolic when sufficiently small symmetry 
breaking terms are added. However, +£2 and —£2 will generically move off 
the X2-axis and will no longer be related to each other by symmetry. Generic 
symmetry breaking will have an analogous effect on +£3 and — £3, and will 
also break the circular symmetry of P and move it off the plane X2 = £3 = 0. 

Sufficiently small symmetry breaking will not change the dimensions of the 
stable and unstable manifolds of ±£2, ±£3 and P, but it will destroy the 
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invariant subspaces, and the heteroclinic connections that existed in the sub- 
spaces will either cease to exist or change their nature. We consider the geomet- 
ric effect of symmetry breaking on each of the former heteroclinic connections 
in turn. 

The heteroclinic connections from ±£2 to ±£3 require the coincidence in 
M 4 of one-dimensional and two-dimensional manifolds; these connections will 
be destroyed by a generic symmetry-breaking perturbation. 

The heteroclinic connections from ±£3 to P occur when the two-dimensional 
unstable manifolds of ±£3 intersect the three-dimensional stable manifold of 
P. Depending on the perturbation, we generically expect to see either transver- 
sal intersections between these manifolds (in which case there are, for example, 
at least two robust heteroclinic connections from +£3 to P) or no intersec- 
tions of the manifolds. The special case where the manifolds are tangent can 
also occur in a codimension-one way. In the case of transversal intersections 
of manifolds, we might expect to see heteroclinic tangles and the associated 
complicated dynamics, depending on whether the dynamics elsewhere in the 
phase space permits reinjection of trajectories into the neighbourhood of the 
transversal intersections. 

The heteroclinic connections from P to ±£2 occur when the two-dimensional 
unstable manifold of P intersects the three-dimensional stable manifolds of 
±£2. There is a clear analogy with the case of connections from ±£3 to P 
and the comments about that case apply equally here. 

While small symmetry-breaking terms generically destroy the heteroclinic 
cycle, there will still be an attractor lying close to the original heteroclinic 
cycle (Melbourne [17] shows this in a closely related case). We show below 
that the form of this attractor (e.g., periodic, quasiperiodic, chaotic) depends 
on the nature of the symmetry-breaking perturbations included. In the fully 
symmetric case, the invariant subspaces defined by X2 = and by X3 = 
restrict each trajectory to one quarter of the phase space, but once the re- 
flection symmetries are broken, a single trajectory may explore more of the 
phase space. We are interested in determining the circumstances under which 
trajectories exhibit switching, i.e., make passages near two or more quarters 
of the original heteroclinic attractor. 
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3 Construction of return maps 

We construct and analyse a return map that approximates the dynamics near 
the cycle. The idea is to define local coordinates and cross-sections near ±£2, 
±£3 and P, then determine local maps valid in a neighbourhood of each of 
±£2, ±-£3 and P, and global maps valid in a neighbourhood of each hete- 
roclinic connection. Composing the local and global maps yields the desired 
return map. Different forms for the return map are obtained depending on 
which of the symmetries are broken. In this section we list the different cases, 
but details of map construction are left to the Appendix. The techniques used 
are, for the most part, standard, although modifications are required to allow 
for the existence of continua of heteroclinic connections along some parts of 
the cycle in the fully symmetric case. 

Throughout, we use a small parameter h to control the size of the local 
neighbourhoods (h -C 1), and small parameters e±, e 2 , £3 to control the extent 
to which the symmetries n±, k 2 , K3 are broken. It turns out to be convenient 
to define the return map on a cross-section near +£3. Using local coordinates 
(n, 61, x 2 , £3) near +£3, where coordinates are chosen so that +£3 is at the 
origin and so that the eigenvectors of the linearised flow align with the coordi- 
nate axes in the manner described in the Appendix, we define a cross-section 

H™ = {(r 1 ,6 1 ,x 2 ,&):0<r 1 <h,\x 2 \ = h,\&\<h} 

and then compute the return map, R : iT 3 n — > H™. The same cross-section 
works equally well near —£3 and the maps R we compute in fact approximate 
the dynamics near any of the four possible paths from ±£3 to ±£3. See the 
Appendix for details. 

Since we are interested in trajectories that switch between positive and 
negative values of x 2 and X3, we introduce the notation ±2 and ±3 to indicate 
whether a trajectory visits +£2 or —£2, and +£3 or —£3. In particular, the 
trajectory starts at one of four possible sections specified by iJ 3 n , and we use 
±2 to specify whether x 2 = +h or x 2 = —h (implying that the trajectory 
recently visited +£2 or —E 2 ). We use ±3 to specify whether the trajectory 
is close to +£3 or —£3. When the trajectory next returns to H™, we will be 
interested in whether it visited +£2 or —£2 en route, and whether it returns 
to +£3 or —£3. 
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First, in the case with full symmetry (ei = €2 = £3 = 0), we have: 
R{r u l5 x 2 = ± 2 h, £3) = (fi = ArJ, 0i = 0i + $ - Qlnn, 

X2 = X2,h = # 2 ), (1) 

where A > and $ are constants, 6 = 616263, Q = (e±e2 + e 2 cs + c^ci) / e^e-^, 
and the constants 6{, e$ and c$ are defined in the Appendix. If X3 > initially, 
the trajectory returns to +E3 after visiting ±2-^2; if X3 < initially, the 
trajectory returns to —E3. 

Second, breaking the K2 and K3 symmetries (ei = 0, e 2 7^ 0, £3 7^ 0) we have: 



R(r 1 ,9 1 ,x 2 = ± 2 /i, £3) = (fi = A ? e 3 ±3 -4i|e 2 ±2 A 3 rf 



01 + $1 + $2 + $3 inn 

e 3 

-—In I e 2 ±2 4»rf 8 | 
ei 



iln 

e2 



e 3 ± 3 Ai|e 2 ± 2 A 3 rf 3 |' 51 , 
x 2 = sgn (e 2 ±2 4»rJ 8 ) /i, | 3 = £2) • (2) 



The trajectory visits +E2 or — £^2 en route according to the sign of €2 ±2 ^43f*i 3 , 
and it returns to +E3 or —E3 according to the sign of €3 ±3 Ai | e 2 ±2 A^r] 3 | 1 • 
Third, if we break the ki symmetry but preserve n 2 and K3 (ei 7^ 0, 62 = 
€3 = 0) we have: 



R(ri,0i,x 2 = ±2/1,6) = [xi = eia T + ^ 2 |^3| 2 cos0i, 

yi = eict; + A 2 |f 3 | 52 sin 0i, 
x 2 = sgn(f 2 )/t, £3 = #2) , (3) 
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where 

x 2 = ±2 (a 3 + ei/ 3 (V - ^ In rf 3 , 

x 3 = ±3 ^i + ei/i + $ 3 - ^-lnn - ^-ln|x 2 |^ |x 2 | 51 , 

#i = #i + 3>i + $2 + $3 - — Inn — — In 1x2 1 — — In |x 3 | . 

e3 ei e 2 

The trajectory visits +£3 or —£2 en route according to the sign of £2, and it 
returns to +E3 or —£"3 according to the sign of X3. In this case, these signs are 
the same as the signs of X2 and X3. The map (3) can be simplified by assuming 
that A3 and A\ are order one and dropping the terms proportional to ei in 
the expressions for £2 and X3. This results in an approximate map: 

R(r 1 ,9 1 ,x 2 = ±2/1, &) = (xi = eia r + Ar[ cos (0 1 + $ - Qlnn) , 

yi = eiai + Arisin(^i + $ - Qlnn) , 

»2 =X 2 ,| 3 = ^2), (4) 

where 5 and Q were defined above, and A and <I> are constants as in equa- 
tion (1). 

Finally, when all symmetries are broken the return map is similar to the 
map (3) above, though with definitions of X2 and x 3 that include terms pro- 
portional to €2 and e 3 : 

x 2 = ±2 U 3 + ei/ 3 (V - ^ Inn) ) r[ 3 

+ e 2 ^1 + ei5 3 ^6»i - ^ In n 

x 3 = ± 3 ^i+ei/i^i + $ 3 -^lnn-^ln|x 2 | ) ) |.r 2 |"' 

+ e 3 ( 1 + eigi \ 0\ + $ 3 — — Inn - — hi |x 2 
V V e 3 ei 

#1 = #1 + $1 + $2 + $3 — Inn - — hi |x 2 | — In I £3 1 . (5) 

e 3 ei e 2 



May 7, 2008 



13:39 Dynamical Systems An International Journal kr 



12 Kirk and Rucklidge 

The trajectory visits +E2 or —E2 en route according to the sign of $2, and it 
returns to +£3 or —E3 according to the sign of £3. It might seem that terms 
proportional to ei in £2 and £3 could be dropped, as they were above. However, 
the terms ±2^3?"i 3 and €2 could nearly cancel and likewise ±3v4i |x2 1^ 1 and £3, 
so we do not drop the ei terms. In fact, it turns out that retaining the e\ terms 
is essential for understanding the switching mechanisms. 

It is possible to write down equivalent maps from Hf 1 — > Hf 1 and H 2 n — * H™. 
Note that the radial coordinates (as defined in the Appendix) play no role in 
the return maps, at the order to which we are working. 

4 Analysis of return maps 

Behaviour in the case without symmetry breaking is well understood and sim- 
ple: whenever 5 > 1 and r is small, iteration of map (1) results in progressively 
smaller values of r and so there is an asymptotically stable heteroclinic cycle. 
The signs of X2 and £3 cannot change, owing to the presence of invariant sub- 
spaces, so each trajectory is confined to one quarter of the phase space. For 
the remainder of this section, we will assume 5 > 1. 

4.1 Global bifurcations 

Global bifurcations are a key ingredient for understanding the dynamics of 
the non-symmetric system. In this section, we describe the global bifurcations 
that are most important for our analysis. 

4.1.1 Homoclinic bifurcation of P. The periodic orbit P has stable and un- 
stable manifolds of dimension three and two, respectively, meaning that trans- 
verse intersections of the manifolds, when they occur, do so in a codimension- 
zero way, while tangencies between the manifolds will be of codimension one. 
Transverse homoclinic orbits can only occur when all symmetries are broken, 
as the following argument shows. If €2 = 0, the subspace %2 = is invari- 
ant; since W S (P) lies in that subspace it cannot intersect W U (P). Similarly, 
if £3 = 0, the subspace X3 = is invariant; since W n (P) lies in that subspace 
it cannot intersect W S (P). If e± = then the rotation symmetry ensures that 
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any intersection of W U (P) and W S (P) will not be transverse. 

In the case ei = 0, e 2 7^ 0, €3 7^ 0, non-transversal homoclinic orbits of 
P occur when one branch of the stable manifold of P is coincident with one 
branch of the unstable manifold of P. This event can be located by calculating 
the image of W U (P) under ^io^ovI/^o^oiI/-^ (see Appendix for definitions 
of the maps fa and Vl/jj) and setting the x 2 component of the image to zero; we 
find that for small symmetry-breaking, non-transversal homoclinic bifurcations 
of P occur at 

e 2 = -± 2 A 3 A 5 J\e 3 \ 5 ^, e x = 0. (6) 

Homoclinic orbits can be formed by coincidence of either of the two branches 
of W U (.P) with either of the two branches of W S (P), resulting in four possible 
homoclinic bifurcations corresponding to the four separate curves implicit in 
the expression above. These curves are shown as dashed lines in figure 3. The 
homoclinic orbit corresponding to the curve in the second quadrant of the 
(€2,63) plane arises from the choice ±2 = + and 63 > 0, and passes close to 
+E 2 and +-E3; the three other bifurcation curves correspond to homoclinic 
orbits with the three other routes past the equilibria, in the obvious way. 

As €\ changes from zero, each curve of non-transversal homoclinic bifur- 
cations will generically split into two curves of homoclinic tangencies, with 
the region between the tangencies being parameter values for which there are 
transverse homoclinic orbits of P. Four curves of homoclinic tangencies and 
two regions of homoclinic tangles are shown schematically in figure 4. 

Inspection of the expression for the x 2 component of the image of W n (P) 
under ^3io</) 3 o\I'23 </ ) 2 oX I'i2 gives more information about loci of the homoclinic 
bifurcations of P when e\ 7^ 0. This component can be written as: 

x 2 = ± 2 i?i 3 (A 3 + ei/ 3 (0i)) + e 2 (1 + £153(61)) , (7) 

where R\ and ©i are complicated functions of the coefficients and parameters. 
In this expression, ^3 + £1/3(81) must remain positive, as explained in the 
Appendix, and R\ is positive. Expressions for the positions of the homoclinic 
tangencies in parameter space can be calculated by setting x 2 = 0; these 
expressions are not included here due to their extreme ugliness. Nonetheless, 
we note that for e\ small, when ±2 = +, there are only solutions with e 2 < 0; 



May 7, 2008 13:39 Dynam 



14 




Figure 3. Schematic bifurcation set for the case ei = 0, ei and 63 small. Regions of asymptotically 
stable quasipcriodic solutions arc bounded by codimension-one curves of global bifurcations, i.e., 
non-transverse homoclinic bifurcations of P (dashed curves) and hctcroclinic bifurcations of the 
cycles ±E2 — > ±£3 — > ±£2 (solid curves). The shapes of the global bifurcation curves correspond 
to the choice 5i > 1, 82 > 1 and S3 > 1, but similar figures could be drawn for the other cases. As 
explained in Section 4.2, the various shading styles indicate the regions in which four different 
quasipcriodic solutions occur. Close to the €2 and £3 axes, two different quasipcriodic solutions 

coexist. 



this is consistent with figures 3 and 4, in which each bifurcation curve is 
confined to a single quadrant. However, if e± is large enough that 1 + 6133(61) 
can change sign as 61 varies, the loci of the homoclinic bifurcations of P can 
change quadrants. Of course, this effect is outside the range of validity of the 
return maps we have constructed, but the principle is worth bearing in mind 
as it appears to influence the dynamics observed in the numerical example 
discussed in Section 5. 



4.1.2 Heteroclinic bifurcation ±E 2 — ► ±E 3 — >■ ±E 2 . In the case that 
all symmetries are broken, consideration of the dimensions of the stable and 
unstable manifolds of the equilibrium points shows that the heteroclinic cy- 
cle +E2 — > +£3 — > +E2 will occur in a co dimension- two manner. However, if 
e\ = 0, the connection +E2 — > +E3 is robust and the intersection of W u (+E^) 
and W s (+£^2) is a codimension-one phenomenon, meaning that the hetero- 
clinic cycle as a whole occurs with codimension one. This latter case is of 
interest since, as we will see, the heteroclinic bifurcation unfolds when ei 7^ 
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Figure 4. Schematic diagram showing part of the bifurcation set for the case ei fixed and non-zero 
but small (compare with figure 3). Dashed curves correspond to homoclinic tangencies of P, solid 
curves in the first (resp. fourth) quadrant correspond to heteroclinic tangencies between W u (— E3) 
(resp. W u (+E3)) and W s (+i?2); and the shading shows regions in which the corresponding 
homoclinic or heteroclinic tangles exist. The dotted horizontal line indicates a path through 
parameter space discussed in Section 4.4; the labels A - D indicate schematically parameter values 



into homoclinic bifurcations of +E2 and +E3 and heteroclinic tangencies be- 
tween yV u (+Ez) and W s (+£ , 2) similar to the way each non-transverse homo- 
clinic bifurcation of P splits into two homoclinic tangencies when ei is varied 
from zero (see above). An analogous argument works for heteroclinic cycles 
involving — E2 and/or — E3. 

Calculations with the local and global maps yields an expression for the 
parameter values at which these heteroclinic bifurcations occur: 



See figure 3. This expression is valid for all four cycles ±£2 — ► ±^3 — > ±£b 
so long as ±3 and the sign of £2 are chosen appropriately. 

4.1.3 Homoclinic bifurcations of ±E 2 and ±E 3 . The dimensions of the 
stable and unstable manifolds of 1L-E2 and are such that if homoclinic 
bifurcations of these equilibria occur, they are of codimension one. 
An argument similar to that used in subsection 4.1.1 shows that we require 



used in section 5. 



€ 3 = - ±3 Al| €2 



ei = 0. 
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€\ 7^ and £3 7^ if a homoclinic bifurcation of ±£2 is to occur, although e 2 
could be zero. Similarly, existence of a homoclinic bifurcation of ±£3 requires 
ei 7^ and e 2 7^ 0, although 63 could be zero. The homoclinic bifurcations 
of ±E2 (resp. ±£3) will be of Shil'nikov type if 82 < 1 (resp. S3 < 1) and if 
c 2 < 2 (resp. c 3 < 2) [34]. 

We can in principle calculate parameter values at which these homoclinic 
bifurcations occur, but the expressions are too nasty to be useful. Instead, we 
note that there can be two homoclinic bifurcations of +E2, one for each branch 
of the unstable manifold of +E2, and a further two homoclinic bifurcations of 
— E2. Similarly, there can be two homoclinic bifurcations of +E% and two ho- 
moclinic bifurcations of —£3. These eight homoclinic bifurcations will in gen- 
eral occur at different parameter values, but in the limit ei — ► 0, will converge 
pairwise on the loci of the four heteroclinic bifurcations involving ±E2 and 
±E% discussed in the previous subsection. For instance, as e\ — > 0, a homoclinic 
orbit of +E2 passing near —E3 and a homoclinic orbit of —E% passing near +£ 2 
will converge in phase space on the heteroclinic cycle +E2 — > — E$ — ► +E2, and 
the parameter values at which the homoclinic bifurcations occur will converge 
in parameter space on the locus of the heteroclinic bifurcation. For clarity, 
these bifurcation curves are not shown in figure 4. 

The dynamics associated with these bifurcations will be discussed further 
below. 

4.2 Breaking the two reflection symmetries 

Here we show that, for e± = and for sufficiently small e 2 and €3, map (2) 
generically has at least one asymptotically stable closed invariant curve and 
the corresponding flow has quasiperiodic solutions. This is not a surprising re- 
sult, since the coordinate 6\ decouples from the other coordinates when e\ = 0, 
in which case our system can be reduced to a three-dimensional system with 
a SSHC between equilibria; earlier work on a system related to our reduced 
system showed that breaking the reflection symmetries can give rise to asymp- 
totically stable periodic solutions [19]. Our main aim in this section is to locate 
the regions in parameter space in which the quasiperiodic solutions exist, for 
comparison with the location of some of the global bifurcations described in 
section 4.1. 
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Figure 5. Schematic graphs of F(ri) (equation (8)) for the choice ±2 = ±3 = 1 when 81 > 1, 
82 > 1, 83 > 1, and A\ = A2 = A3 = 1. The solid (resp. dotted) curve indicates values of r\ for 
which the next values of ±2 and ±3 are (resp. are not) both positive; we seek values of n for which 
the solid curve intersects the diagonal. For small positive £2 and £3, a stable fixed point exists (see 
panel (b)). This fixed point ceases to exist in the second quadrant when e2 = — A3A2 2 eg 2 * 53 , when 
there is a non-transversal homoclinic bifurcation of P (limiting case shown in panel (a)). The fixed 
point is destroyed in the fourth quadrant when ez = — Aie* 1 , when there is a non-transversal 
heteroclinic connection from +£3 to +E2 (limiting case shown in panel (c)). 



The r± component of map (2) is independent of the other variables, and so 
we first seek values of r± for which F(r±) = n, where 



F( ri ) 



A, 



8 3 \°i 



(8) 



For each choice of £2 and £3, there are two possible signs of each of ±2 and ±3, 
but the case (±2 = +, £2 > 0) is equivalent to (±2 = — , £2 < 0), and the case 
(±3 = +, £3 > 0) is equivalent to (±3 = — , £3 < 0). Without loss of generality, 
we focus on the case ±2 = +, ±3 = +, and seek values of £2 and £3 for which 
there exist fixed points of map (8). Fixed points of this type have positive 
values within the absolute value signs in (8), since the signs of £2 + A^rf 3 and 
£3 + Ai|£2 + ^43r^ 3 |' 51 determine the next values of ±2 and ±3. 

For sufficiently small, positive £ 2 and £3, F(0) = ^(£3 + Aie^ 1 ) 52 > 0. For 
r\ larger than £2 and £3 but still smaller than one, we have F{?\) ~ rf , which 
is less than r\ since 5 > 1. Thus, by continuity, the map has a fixed point (see 
figure 5(b)). Since F(r±) is monotonically increasing, the slope of F at the 
fixed point is less than one, so a stable fixed point exists for £2 > 0, £3 > 0. 

This fixed point (i.e., a fixed point with ±2 = ±3 = +) also exists in parts of 
the second and fourth quadrants of the (£2, £3) parameter plane. To determine 
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the region of existence in the fourth quadrant, we fix £2 at some small positive 
value and decrease €3. This shifts the graph of F(r±) down, from which it is 
found that a stable fixed point exists until ^(0) = 0, i.e., until €3 = —Aie^ 1 
(figure 5(c)). Thus the fixed point ceases to exist in the fourth quadrant at 
the locus of the heteroclinic bifurcation from +E3 to +E2 (c.f. section 4.1.2). 

To determine where the fixed point exists in the second quadrant, we fix 63 at 
some small positive value and decrease €2- This decreases F(0) and also changes 
the shape of the graph of -F(n); the graph remains monotonic increasing in 
r\ while €2 is positive, but develops a turning point once €2 becomes negative, 
with F(r\) decreasing for r\ near zero. The decreasing section is indicated by 
a dotted curve in figure 5(a), and corresponds to future values of ±2 and ±3 
not both being positive. For small enough negative €2 the dotted section of the 
graph lies to the left of fixed point, but when €2 = — AzA 2 2 e 3 2 3 , the dotted 
curve reaches to the diagonal and the fixed point ceases to exist (figure 5(a)). 
Thus the fixed point ceases to exist in the second quadrant on the locus of the 
homoclinic bifurcation of P (c.f. section 4.1.1). 

The region of existence of this stable fixed point is indicated by the left- 
leaning close hatching in figure 3. Calculations with the other combinations of 
signs of ±2 and ±3 are analogous, and yield different regions of existence for 
the corresponding fixed points. Fixed points may coexist as shown in figure 3. 

Since for fixed r\ the 6\ component of map (2) is a rigid rotation, a fixed 
point of equation (8) generically corresponds to a closed invariant curve in 
(2) and to a quasiperiodic solution in the full flow. The angle 9\ decouples 
from the rest of the dynamics, and so the full flow will have an invariant torus 
foliated by periodic orbits for a dense set of parameter values. Stability of 
these solutions follows from the stability of the fixed point of (8). 

The calculations above were for the case 8\ > 1, 82 > 1, S3 > 1. Similar 
calculations done when one or more of the <5j is smaller than one lead to 
similar regions of existence of quasiperiodic solutions, except that there are 
additional saddle-node bifurcations of the tori close to the relevant global 
bifurcations; these saddle-node bifurcations arise since the global bifurcations 
destroy solutions of different stabilities depending on the sign of Si — 1. 

The special case that precisely one of €2 and €3 is zero (i.e., only one reflection 
symmetry is broken) is covered by the analysis above; there will be two fixed 
points of the map with corresponding (foliated) tori in the flow. 
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In summary, when the two reflection symmetries are broken but the ro- 
tation symmetry is preserved, the flow generically has asymptotically stable 
quasiperiodic solutions that do not exhibit switching. Different quasiperiodic 
solutions coexist in regions bounded by global bifurcations that will play an 
important role in generating complex dynamics when e\ ^ 0. 

4.3 Breaking the rotation symmetry 

In this subsection we show that for sufficiently small values of e±, with €2 = 
£3 = 0, the map (4) has a stable fixed point. 

By rescaling r\ and e\ by order one amounts and moving the origin of the 9\ 
coordinate, we can without loss of generality set A = 1, a r = 1 and a; = in 
(4). Ignoring for now the X2 and £3 components of the map and working with 
polar coordinates (p,4>) centred at (xi,y±) = (ei,0) (so that x\ = e± + pcos(f> 
and y\ = psmfi), map (4) reduces to 

X\ = e± + /ocos (p 
< _ where 

yi= psm(j) 

The constant <3? may take a different value here than in equations (4). Fixed 
points of (9) satisfy r\ = \Jx\ + yf and 6\ = arctan (yi/xi). To find solutions 
of the first of these equations, note that circles of radius r about (xi,yi) = 
(0,0) map under (9) to circles of radius r s about (ei,0). Since 5 > 1 and for 
small r, these circles will intersect if r > e\—r 5 and r < ei+r 5 ; the intersection 
points are candidate fixed points of the map. For each small fixed value of e±, 
there will be some non-zero interval a < r < b on which the inequalities are 
both satisfied. See figure 6(a-c). 

The second equation, 6\ = arctan (yi/xi), is satisfied for at least one value 
of r in [a, b], as the following argument shows. When r = a, the circles 
(j"i, #i) = (r, 61) and (p, 0) = (r 6 , 0) intersect at a single point, (ri,9\) = (a, 0), 
alternatively (p,</>) = (a s ,ir). As r is increased beyond a, the intersection 
point splits into two (with corresponding (f> values just below it and just 
above — it). The intersection points come together again at (r±,6i) = (6,0) or 
= (b s ,0), in the manner shown in figure 7(a). The corresponding value 
of <f> for each intersection point can be calculated from (9) (see figure 7(b)) 



= 6»i + $ - Qlnn. 



(9) 
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(a) (b) (c) 




Figure 6. Finding fixed points of equations (9). Panels (a-c) show schematically the relative 
positions of the circle of radius r around (xi, yi) = (0, 0) (large circle in each panel) and its image 
under map (9) for various sizes of r: (a) r < ei — r*, no intersection; (b) r > e\ — r s and 
r < e\ + r s , one or two intersections; (c) r > e\ + r s , no intersections. 




b r a b r 



Figure 7. (a) Schematic <j> values for intersections of the two circles (ri,Q\) = (r, 8\) and 
(p,<t>) = (t* 5 , , in the situation shown in figure 6(b), plotted as a function of r; (b) Schematic 
showing 8\ (solid and dashed curves) and — Qlnri (dotted curve) for points of intersection of the 
two circles, plotted as a function of r at the intersection points; (c) <j> and <j> values at points of 
intersection of the two circles. Solid curves correspond to the upper intersection points and their 
images, dashed curves correspond to the lower intersection points and their images. 

from which it is seen that the two branches of (j), arising from the upper and 
lower intersections of the two circles, start and end at the same point as each 
other, as shown in figure 7(c). At least one of the two branches of the graph 
of 4> vs r therefore intersects the graph of 4> vs r for that same branch for at 
least one r in [a, b]. Thus, there is at least one fixed point of the map (9). 

It is straightforward to show that the determinant of the Jacobian of map (9) 
is <5r 2 ( 5-1 ) and that the absolute value of the trace of the Jacobian is bounded 
above by \/ (S + l) 2 + Q 2 r 5_1 . For each choice of e±, the value of r at the 
corresponding fixed point lies in [a, b], where a and b depend on e\ with a > 
and a and b tending to zero as ei — > 0. Thus, for sufficiently small ei, the 
determinant and trace of the linearised map are small enough to ensure that 
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the relevant fixed point is stable. It follows that the corresponding fixed point 
of equations (4) is also stable, since orbits of that map collapse onto a constant 
value of the third (radial) coordinate after one iteration of the map. 

In summary, for sufficiently small ei when 62 = £3 = 0, the heteroclinic cycle 
that occurred in the fully symmetric case is replaced by a stable periodic orbit. 
Using the reflection symmetries K2 and K3, we find four stable periodic orbits 
co-exist, one in each quarter of the phase space. 

4.4 Breaking all symmetries 

Direct analysis of the return map valid for the case that all symmetries are 
broken is not feasible because of the extremely complicated form of that map. 
Instead, in this section we use our knowledge of the dynamics in the case 
ei = and arguments about generic unfoldings of this special case to deduce 
what types of dynamics will be seen in the fully asymmetric case for ei near 
zero, and to (approximately) locate each type of behaviour in parameter space. 
This procedure allows us to make specific predictions about the mechanisms 
underlying the complicated dynamics observed in numerical examples, such 
as the example described in section 5. We are especially interested in finding 
mechanisms that cause repeated, non-periodic switching in our system. 

The dynamics associated with the case ei = is summarised in figure 3, 
which shows eight curves of global bifurcations bounding regions in which 
there are asymptotically stable quasiperiodic solutions. In this case the rota- 
tion symmetry prevents coupling of the two frequencies associated with each 
quasiperiodic solution and the dynamics is simple. Once ei moves away from 
zero, we will generically see locking of the frequencies. For instance, if we were 
to fix ei sufficiently small but non-zero and pick €2 and £3 positive and with 
values midway between the two pairs of global bifurcation curves in the first 
quadrant of figure 4, then along a one-dimensional path through the parameter 
space such as the dotted line in figure 4 there will be intervals of quasiperiod- 
icity interspersed with intervals of locked, periodic behaviour. Associated with 
the frequency locking there may be complicated dynamics such as period- 
doubling cascades and chaotic dynamics. However, this behaviour will mostly 
be confined to regions of phase space near the original quasiperiodic solutions, 
and is not the main mechanism for switching in our system. 



May 7, 2008 



13:39 Dynamical Systems An International Journal kr 



22 Kirk and Rucklidge 

As shown in section 4.1.1 and in figure 4, each curve of non-transverse ho- 
moclinic bifurcations of P seen in figure 3 turns into a wedge in parameter 
space of transverse homoclinic orbits of P when ei changes from zero. There 
will be horseshoes and chaotic dynamics associated with the transverse homo- 
clinic orbits, although the chaos may not be attracting. In numerical examples 
we might expect to see a mixture of stable periodic orbits and stable chaotic 
dynamics, in overlapping regions of parameter space. 

An interesting consequence of the occurrence of homoclinic tangles is that it 
provides a mechanism for switching of orbits with respect to the X2 variable. 
For instance, for sufficiently small ei, and with €2 > 0, £3 > and both 
small, solutions that make excursions near +E2 can get trapped. The trapping 
region is bounded in part by one branch of W S (P) and trapped solutions make 
excursions past +E2 but cannot cross W S (P) so cannot get close to —E2. If 
62 is decreased, say by moving along the dotted path shown in Figure 4 into 
the homoclinic wedge in the second quadrant, the trapping region develops a 
leak when a homoclinic tangency forms between that branch of W S (P) and 
a branch of W U (P); solutions are then able to cross W S (P), and may visit a 
neighbourhood of —E2. We call this 'switching in X2 '■ Switching of this type 
(from positive to negative X2) can occur for parameter values to the left of the 
right boundary of the homoclinic wedge in the second quadrant. A numerical 
example of this leaking process is given below in figure 9. 

Once a switching solution arrives in the region with X2 < 0, it may then 
get captured by an attractor lying solely in the negative X2 region of phase 
space, in which case no more switching will be observed. Alternatively, if there 
is a mechanism for orbits to leak back to the original region of phase space 
then there could be sustained switching in X2- This latter case cannot occur 
for arbitrarily small e\ as the following argument shows. For the case ei = 0, 
results from section 4.2 show that orbits which make repeated excursions near 
—E2 and +-E3 occur in the second quadrant of figure 4 and in the first quadrant 
as far as the locus of homoclinic bifurcations of P. (This homoclinic bifurcation 
involves a different branch of W U (P) than the homoclinic bifurcation occurring 
for €2 < discussed in the last paragraph.) For small ei, there will be a wedge 
of homoclinic tangencies of the relevant branches of W n (P) and W S (P), with 
the wedge lying entirely within the region £2 > 0, e% > 0; orbits that make 
a number of excursions past — E2 before switching and passing close to +E2 
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can only occur for values of £2 an d £3 lying to the right of the left boundary 
of the wedge. Thus, for sufficiently small e±, there is no overlap between the 
region of parameter space where there is switching from positive to negative 
X2 and the region where there is switching from negative to positive X2, with 
the consequence that there can be no sustained switching in x%- 

However, as argued in section 4.1.1, for large enough t\ the curves of ho- 
moclinic tangency may change quadrants in the (€2,63) parameter plane, and 
then the switching regions can overlap, making sustained switching in %2 pos- 
sible. As pointed out in section 4.1.1, our return map construction is not valid 
for 'large' e±, so we have not proved the existence of sustained switching, just 
shown how it might feasibly occur. It is not possible to determine a priori 
how big ei would have to be to get sustained switching, but we have shown 
that there is a threshold in ei below which sustained switching in %2 is not 
possible. The value of this threshold does not go to zero as £2 and £3 go to 
zero, as it comes from the requirement that 1 + £153(61) changes sign (as a 
function of Oi) in equation (7). Another way of understanding this is to note 
that in (7), if ±2 = + and £2 > 0, then the only way of having X2 negative 
is to have 1 + £153(61) < for some value of ©i. This is a necessary but not 
sufficient condition, as the attractor may not explore the required range of 61. 

The four curves of heteroclinic bifurcations of the cycles ±£2 — ► ±£3 — > 
±£2 shown in figure 3 also split when £1 becomes non-zero, being replaced by 
eight curves of homoclinic bifurcations and eight curves of heteroclinic tangen- 
cies between W U (±E^) and W S (±E2), as described in section 4.1.2. If they are 
of Shil'nikov type, the homoclinic bifurcations can complicate the dynamics by 
inducing chaotic dynamics. The heteroclinic bifurcations are associated with 
switching in the 23 coordinate similarly to the way switching in X2 is associ- 
ated with homoclinic bifurcations of P, described above. More precisely, the 
eight curves of heteroclinic tangencies between W u (±Es) and W s (±i?2) come 
in pairs, with each pair bounding a wedge in parameter space. At parameter 
values within each wedge there is a heteroclinic tangle of one pair of manifold 
branches. For instance, for sufficiently small £1 there will be a heteroclinic 
wedge involving one branch of W u (+£ , 3) and one branch of W s (+E2) occur- 
ring in the fourth quadrant in the (£2,63) plane. Above this wedge, W s {+E2) 
bounds in part a trapping region; orbits in the trapping region make excur- 
sions past +£3 but cannot cross W s (+E2) to get close to — £3. The trapping 
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region develops a leak when a heteroclinic tangency forms between the appro- 
priate branches of W u (+£3) and W s (+£ , 2) thus allowing solutions to cross 
W s (+E2). We call this 'switching in X3'. An argument analogous to that used 
for switching in X2 can be used here to show that there is a (generically dif- 
ferent) threshold in ei below which there can be no persistent switching in X3. 
This threshold does not go to zero when €2 and £3 go to zero. 

The mechanisms inducing switching in X2 and in X3 are distinct, but or- 
bits that switch persistently in both x<i and X3 are possible for ei above the 
thresholds for both mechanisms. Switching in each variable requires the rota- 
tion and appropriate reflection symmetry to be broken. It is possible to have 
persistent switching in X2 with 63 = 0, or switching in X3 with 62 = 0, though 
we will not explore this possibility in detail. The example in section 5 shows 
that persistent switching is easily observed in numerical examples. 

Both of the global bifurcations we have identified as inducing switching, i.e., 
homoclinic tangencies of W U (P) and W S (P) and heteroclinic tangencies of 
W u (±£^3) and W S (±E2), will produce horseshoes in the dynamics. In the case 
of homoclinic tangencies, this is a standard result and in the case of the hetero- 
clinic tangencies, reinjection into the neighbourhood of the heteroclinic tangle 
is provided by proximity in phase and parameter space to the heteroclinic cy- 
cle ±£2 — * ±£3 —* ±-#2- I n either case, we expect the onset of switching to 
be commonly associated with nearby chaotic dynamics; chaotic orbits before 
the onset of switching, chaotic transients for switching orbits and orbits that 
switch chaotically might all be seen. However, other types of switching are 
also possible, such as periodic switching where the attractor is a periodic orbit 
that crosses the (non-invariant) hyperplanes X2 = and X3 = or periodic 
switching where the attractor is a 'noisy periodic orbit' such as results from a 
cascade of period doubling. In the latter case the itinerary of visits to ±£2 or 
±£3 will be periodic even though the actual orbits are not. 
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5 Example 

We consider the following system of equations to illustrate the dynamics of 
interest in this paper: 

zi = (1 + i)zi - {z^zi - (c 2 + l)xjzi + (e 3 - l)x\zi + eidn + eid i2 xi , 

(10) 

x 2 = x 2 - x 2 - (c 3 + l)x§x 2 + (ei - l)|zi| 2 x 2 + e 2 d 2 i + ei<i 22 xix 2 + eie 2 d 23 xi , 

(11) 

x 3 = x 3 - xl - (ci + l)|zi| 2 x 3 + (e 2 - l)x2X 3 + e 3 cZ 3 i + ei<i 32 xix 3 + eie 3 <i 33 xi , 

(12) 

These equations were derived by starting with the structurally stable hete- 
roclinic cycle considered in [2] , turning a pair of equilibria of that cycle into 
a periodic orbit by adding a trivial phase variable, and adding the simplest 
possible terms that break the symmetries in generic ways. The parameters ei, 
e 2 and e 3 in equations (10-12) play the same role as in the maps derived earlier 
in this paper. 

The model used in [20] is similar, but differs in two respects. First, the 
symmetry-breaking terms in [20] were fifth-order in the x and z variables, 
rather than constant, linear and quadratic here. Second, the model in [20] 
respects the symmetry (zi,x 2 ,x 3 ) — > (— z±, — x 2 , — x 3 ), as appropriate for a 
model of a dynamo instability: the invariant subspace z\ = x 2 = x 3 = 
corresponds to the absence of any magnetic field. We do not expect the first 
difference between models to alter the qualitative behaviour, but the enforced 
symmetry may have a significant effect, as discussed briefly in the next section. 

The coefficients in the equations were chosen to be: ci = 1.2, e\ = 1.0, c 2 = 
1.1, e 2 = 1.0, c 3 = 1.1, e 3 = 1.0 for the contracting and expanding eigenvalues, 
and du = du = 10- 4 , d 2 i = 10" 1 , d 22 = 10" 1 , d 23 = 10 3 , d 31 = 10~ 3 , 
^32 = 10 -4 , d 33 = 1 for the symmetry-breaking coefficients. The eigenvalues 
were chosen to be of order one, with contraction dominating expansion at each 
point (Si = 1.2, 5 2 = 1.1, <5 3 = 1.1, and an overall 5 = 1.452). The symmetry 
breaking coefficients are notionally small, but those coefficients (d 23 and d 33 ) 
that are multiplied by two e's were chosen to be larger to compensate for this. 
The exact numbers are not important, though they will affect the details of 
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what is observed. However, choosing c?23 and ^33 to be reasonably large means 
that the switching dynamics is easier to obtain for small values of e±: in order 
to get persistent switching, the €162^23X1 and €163^33X1 terms in (10-12) need 
to be reasonably important. 

We integrated the equations numerically using the Bulirsch-Stoer adaptive 
integrator [35], with a tolerance for the relative error set to 1CT 12 for each 
step. Poincare sections were computed using algorithms from [36]. 

By varying ei, €2 and £3, we are able to find examples of the important 
symmetry-breaking effects discussed in the previous sections of this paper. 
The cases with full or partial symmetry preserved give straightforward results, 
which we describe only briefly; more details are provided of the case of fully 
broken symmetry. 

If all symmetries are preserved (all ej = 0) then each solution starting off the 
invariant subspaces is attracted to one of four symmetry-related structurally 
stable heteroclinic cycles. If ei 7^ 0, 62 = £3 = (rotation symmetry broken, 
reflections preserved), numerics confirm the predictions of section 4.3, and a 
single attracting periodic orbit is found in each quarter of the phase space. 
If rotation symmetry is preserved as well as one reflection, and the other 
reflection is broken, then solutions are attracted to a foliated torus, as discussed 
in section 4.2. In the case that both reflections are broken but the rotation 
symmetry is preserved, numerics confirm the predictions of section 4.2; we 
find that there exist attracting quasiperiodic solutions in regions bounded by 
curves of global bifurcations, as shown schematically in figure 3. Analysis of 
the maps derived earlier allows us to predict scaling of the loci of various global 
bifurcations in the limit of small symmetry breaking. For instance, equation 
(6) tells us that for e\ = and £3 — > 0, homoclinic bifurcations of P associated 
with equations (10-12) occur for €2 = constant x | €3 1 ^s* 53 but the value of 
the constant is not determined by the map analysis. Numerical simulations of 
equations (10-12) confirm the scalings for the various global bifurcations. 

To illustrate the phenomena associated with breaking all symmetries, it is 
helpful to consider the changes in dynamics seen along a one-dimensional path 
such as that shown as the dotted line in figure 4. We first chose a value of ei 
below the thresholds for persistent switching in X2 and X3. For instance, fixing 
ei = 10~ 4 and £3 = 0.001 and allowing 62 to vary, we see the following types 
of dynamics. 
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Figure 8. Onset of switching in X2 for equations (10—12), associated with crossing into the wedge 
of homoclinic bifurcations of P in the second quadrant of (e2 , £3)-space. Parameters ei = 10~ 4 , 

€3 = 0.001 are fixed and ei is decreased: (a) €2 = —5.67 X 10~ 7 gives a chaotic attractor confined to 

the region X2 > 0; (b) €2 = —5.68 X 10~ 7 gives a chaotic transient with X2 > 0, then the sign of X2 
changes and the orbit is attracted to a quasiperiodic solution with X2 < 0. Other coefficients as 

described in text. The chaotic nature of the orbit in (a) and the transient in (b) is not apparent on 
the timescale used to plot the time series. 



Picking £2 = 3 x 10~ 5 yields a point (labelled A) lying to the right of the 
homoclinic wedge and to left of the heteroclinic wedge in the first quadrant 
of figure 4. For these £2 and €3 values but for ei = 0, there exists an attract- 
ing quasiperiodic solution with and X3 both positive. With ei = 10 the 
same type of quasiperiodic solution exists. As €2 is decreased while e\ is fixed 
at 10 we find, as expected, intervals of 62 in which there are quasiperiodic 
attractors interspersed with intervals on which there is locking of the two fre- 
quencies associated with the quasiperiodic solution (periodic orbits). In some 
intervals, period doubling cascades are observed, as is normal near quasiperi- 
odic behaviour. The interchange between locking and quasiperiodic behaviour 
persists until we approach the homoclinic wedge at negative values of €2- 

As this homoclinic wedge is approached, apparently chaotic dynamics is ob- 
served, consistent with the appearance of horseshoes associated with the im- 
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Figure 9. Illustration of how leaking develops as the attractor crosses W S (P) (approximately 
X2 = in this figure), (a) Poincare map for the orbit shown in figure 8(a); (b) Poincare map for the 
orbit in figure 8(b). In each case the Poincare section is i?J n with \xg\ = h = 0.005. In each panel, a 

dot (resp. cross) indicates that the orbit next crosses the Poincare section with X2 > 
(resp. X2 < 0). In (b), the upper collection of dots corresponds to the chaotic transient, which ends 
in a single cross, after which the orbit switches to a quasipcriodic attractor, represented by the 

lower collection of crosses. 



pending homoclinic tangency. Before the first tangency is reached (labelled B 
in figure 4), orbits are trapped in the region with X2 > 0, X3 > (figure 8a). 
If 62 is decreased past the tangency value (labelled C), the attractor crosses 
W S (P) and so a typical orbit will display a chaotic transient with X2 > 0, 
X3 > 0, and then switch to X2 < 0, £3 > 0, after which the orbit is attracted 
to a quasiperiodic solution in that quarter of phase space (figure 8b). Cor- 
responding Poincare maps are shown in figure 9. With negative values of €2, 
once the trajectory has switched to X2 < 0, the behaviour is analogous to 
that observed with €2 > and X2 > 0: quasiperiodic attractors interspersed 
with frequency locking. Note there is no persistent switching in X2, and no 
switching in X3, for these parameter values. 

If 62 is now increased from 62 = 3 x 10~ 5 while ei and €3 are kept fixed as 
before, we approach the wedge of heteroclinic connections from —E% to +E2 
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(following the dotted line in figure 4). Once the left edge of the heteroclinic 
wedge has been crossed, orbits can switch from X3 < to x 3 > (but not from 
X3 > to X3 < 0). For example, at €2 = 1.1 x 10 -4 (labelled D in figure 4), 
we find a chaotic transient with £3 < that has a single switch to a locked 
periodic orbit with x 3 > 0, in a manner similar to the single X2 switch in 
figure 8(b). For larger €2, we find other examples of quasiperiodicity, locked 
periodic orbits, chaos, chaotic transients and single switches from X3 < to 
X3 > 0, consistent with the analysis presented above. We did not find any 
examples of persistent switching. 

The behaviour for negative values of £2 and/or £3 is analogous: single 
switches can be found, but there is no persistent switching for e\ = 10~ 4 . 
However, persistent switching in one or both of X2 and x 3 is observed if we 
increase t\. Figures 10 and 11 show an example of persistent switching in %2 
(but not x 3 ) for ei = 3 x 1(T 4 , e 2 = 2 x 1CT 4 and e 3 = 0.001. The trajectory 
crosses the Poincare section in a curve that appears reasonably smooth at 
the largest scale (figure 11), but the magnified inset shows that the curve has 
structure, and that parts of the curve lie below W S (P), leading to switches 
from X2 > to X2 < 0. For these parameter values, there appears to be no 
attractor with X2 < 0, and after a short transient, the trajectory switches back 
to X2 > 0. Dynamics with a larger value of ei = 0.005 is shown in figures 12 
and 13: here we have persistent switching in X2 and x 3 . 

Our understanding of the mechanism behind persistent switching in X2 or 
x% requires that e\ be large enough that 1 + tig 3 {6) or 1 + e\gi{9) can take on 
positive and negative values, as a function of 9. In order to illustrate this effect, 
we have fitted these functions using the trajectories in figures 10 (ei = 3x 10~ 4 ) 
and 12 (ei = 0.005). We have concentrated on switching in x 3 , so the fact that 
these trajectories are for different values of 62 does not affect our conclusions. 
We took the coordinates (x 3 ,6\) on the section H° nt and x 3 from the next 
intersection with H™, and used these data to fit a map of the form of ^12: 

x 3 = x 3 (Ai + ei/i(0i)) + e 3 (1 + ei0i(0i)) , 

see equation (18). We represented the two functions f\ and g\ by a finite 
Fourier series, and were able to fit the data to within one part in 1000 for all 
points with the smaller value of e±, and to within one part in 100 for all but 
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Figure 11. Poincare section corresponding to the time series shown in figure 10. The inset shows 
an enlargement of the region in the box marked in the main picture. The Poincare section is 
with |oS3 1 = h = 0.01. A dot (resp. cross) indicates that the orbit next crosses the Poincare section 
with X2 > (resp. X2 < 0). The inset shows that crosses (indicating a switch) occur where the 
trajectory lies below W S (P) (approximately X2 = 0). 
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(a) 




Figure 12. Time series showing persistent switching in X2 and £3, for (10-12) with ei = 0.005, 
£2 = 3 X 10~ 5 , €3 = 0.001. Other coefficients as defined in text. 



a handful of the points for the larger value of e±. As expected, the A\ + ei/i 
part of the map remains positive, but 1 + eigi can change sign for the larger 
value of ei, as shown in figure 14. Indeed, the numerical ratio of the amplitudes 
of the two fitted functions is 16.663 while the ratio of the two values of ei is 
16.667. Our understanding requires the change in sign of l+e±gi as a necessary 
condition for persistent switching in X3, which is confirmed by this illustration 
and by our other calculations. A similar transition occurs (at a smaller value 
of ei) at the onset of persistent switching in X2, and in that case, the data can 
be fitted to within 1 part in 10,000 or better. 

We note that we were able to find parameter values associated with persis- 
tent switching in x 2 alone but were unable to get persistent switching in x 3 
without also having switching in x 2 - This is a consequence of the particular 
choice of symmetry-breaking coefficients we use: the threshold in ei for per- 
sistent switching in X3 is higher than the threshold for persistent switching in 
X2 for the chosen coefficients (since we have (I23 > ^33). For other parameter 
choices, the thresholds could be the other way around. 
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Figure 13. Poincare section corresponding to the time series shown in figure 12. The inset shows 
an enlargement of the region in the box marked in the main picture. The Poincare section is _ff} n 
with \x$ \ = h = 0.01. Four symbols are used: a dot (resp. cross) indicates that the orbit next 
crosses the Poincare section with X2 > (resp. X2 < 0) and with X3 > 0. A + (rcsp. square) 
indicates that the orbit next crosses the Poincare section with X2 > (resp. X2 < 0) and with 
X3 < 0. The division between orbits falling either side of W S (P) is clearly visible. Orbits falling on 
opposite sides of W s (±i?2) are reasonably well mixed with this choice of cross-section. 
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Figure 14. The function 1 + eigi(6) fitted to the data in figures 10 (ei = 3 x 10~ 4 ) and 12 
(ei = 0.005). This function must be both positive and negative as a function of 8 in order to allow 

persistent switching in X3. 

In the numerical simulations described above, the quantities 5\, 62, S3 and 
5 were all greater than one. This choice was made to ensure that that the 
heteroclinic cycle was attracting in the fully symmetric case and to remove 
any possible complications due to chaotic dynamics associated with homoclinic 
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bifurcations of =Li?2 an d ±E$. Much of the same switching dynamics will still 
occur if one or both of 82 and £3 is less than one while 5 > 1, since the 
mechanism for switching we have found does not depend on the size of the 
individual 5i. However, in this case there may be additional complications in 
the dynamics associated with the homoclinic bifurcations of the equilibria. 

6 Conclusions 

This paper has investigated the effect of small symmetry-breaking on the dy- 
namics near a structurally stable heteroclinic cycle connecting two equilibria 
and a periodic orbit. The heteroclinic cycle is structurally stable in the case 
that there are two reflection symmetries and a rotation symmetry in the un- 
derlying system; we were interested in the dynamics seen when one or more of 
the symmetries is broken. It was reported in [20] that this type of system can 
exhibit seemingly chaotic dynamics along with repeated but irregular switch- 
ing of sign of various variables, but details of the mechanisms underlying the 
onset of complicated dynamics were not explored there. In this paper, we have 
identified global bifurcations that induce the onset of chaotic dynamics and 
switching near a heteroclinic cycle of this type. These turn out to be homo- 
clinic tangencies between the stable and unstable manifolds of the periodic 
orbit, and specific heteroclinic tangencies between stable and unstable mani- 
folds of the two equilibria. By construction and analysis of approximate return 
maps, we were able to locate (approximately) the global bifurcations in pa- 
rameter space and hence to isolate instances of the different types of switching 
and chaotic dynamics in a specific numerical example. 

In addition to identification of the mechanisms underlying the onset of 
switching, two important insights have been gained from this study. First, we 
found that interaction of the different symmetry-breaking terms is required 
for switching; partial symmetry breaking (where one or two of the three sym- 
metries are retained) did not result in switching. Switching results from the 
right combination of a global bifurcation (which results in turn from breaking 
of the rotation symmetry) and small breaking of at least one of the reflection 
symmetries. Second, we found there is a threshold in ei below which there can 
be single switches in the signs of certain variables but no persistent switching. 
The important point here is that persistent switching does not result from ar- 
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bitrarily small symmetry breaking, but is a 'large' symmetry-breaking effect. 
Of course, 'small' and 'large' are relative terms, and addition of seemingly 
tiny symmetry-breaking effects might actually result in persistent switching, 
as was the case in the numerical example we investigated in section 5. 

One aspect of this problem which has not yet been investigated is whether it 
is possible to make a priori predictions about switching rates or derive scaling 
laws for switching times. It is plausible that switching rates and times might 
depend on the 'distance' from the global bifurcation that induces the switching, 
but no detailed attempts have yet been made to quantify such a relationship. 
The statistics of switching intervals were measured in the related model of [20] , 
who report an exponential distribution of intervals between switches. 

Finally, we note that the dynamo model in [20] has a symmetry that is never 
broken (this is the symmetry (zi,X2,x 3 ) — > (— z±, — X2, — £3) in the notation 
of [20]). Retention of this symmetry while breaking all others amounts to 
retaining invariance of the z\ = X2 = x 3 = subspace, and will have a 
consequence of relating the dynamics in different parts of the phase space. For 
example, if it is possible to switch from (x2 > 0, X3 > 0) to (x2 > 0, X3 < 0), it 
will also be possible to switch from (x2 < 0,X3 < 0) to [x2 < 0,x 3 > 0). Our 
results do not include this effect, and retaining this symmetry may well have 
profound effects on the switching properties. Nevertheless we expect our basic 
ideas about switching being induced by a balance between a global bifurcation 
and symmetry-breaking terms and the existence of a threshold for persistent 
switching to apply quite generally, and to the example in [20] in particular, 
even if the details turn out not to be directly relevant. 



7 Appendix: Details of return map construction 
7.1 Coordinates and cross-sections 

Following [4], we distinguish radial, contracting, and expanding directions near 
the equilibria in the fully symmetric case. If V\ = {(zi,X2,x 3 ) : x 3 = 0}, 
V2 = {(zi,x 2 ,x 3 ) : zi = 0}, Vz = {(zi,x 2 ,x 3 ) : x 2 = 0}, with V = V 3 , 
then the radial eigenvalues at ±Ej (j = 2, 3) are the eigenvalues of the lin- 
earised vector field at ±Ej (i.e., eigenvalues of (df)±E j ) restricted to VjOVj-i- 
The contracting eigenvalues are the remaining eigenvalues of (df)±E j in Vj-i, 
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and the expanding eigenvalues are the remaining eigenvalues in Vj. The ra- 
dial direction is then the span of the eigenvectors corresponding to the radial 
eigenvalues, and similarly for the contracting and expanding directions. Near 
P we define the radial direction to be the direction of V\ n V3 (i.e., the plane 
X2 = £3 = 0), the contracting direction is parallel to the X3-axis, and the 
expanding direction is parallel to the X2-axis. These definitions are consistent 
with those in [4] but are adapted for the case where there is a periodic orbit 
in the hetero clinic cycle. 

We choose local coordinates near each of P, ±i?2j and ±E^ to make the 
linearised dynamics as simple as possible. Near +E2 in the fully symmetric 
case, we define £2 = £2 — X2, where X2 is the value of X2 at +E2, and then 
use local coordinates (zi, £2, £3); zi, £2 and X3 correspond to the contract- 
ing, radial and expanding directions, respectively. Under symmetry breaking, 
+E2 moves in proportion to the magnitude of the symmetry breaking, and 
the local coordinates are measured from the new position of the equilibrium 
point. The eigenvalues and eigenvectors change similarly, but since the eigen- 
values are generically distinct and non-zero, small symmetry-breaking will not 
change the nature of the local structure and we can use the slightly altered 
eigenvectors to define a slightly altered local coordinate system. We continue 
to identify radial, contracting and expanding directions once weak symmetry 
breaking is introduced, in the obvious way, and retain the notation (zi, £2, ^3), 
for the altered coordinates, although z\ and X3 may no longer coincide with 
the corresponding global coordinates. 

A similar construction is used near —E2 except that £2 = —X2+X2, where X2 
is the value of X2 at —E2. The point of defining £2 in this way is that positive 
values of £2 near +E2 are mapped under the reflection K2 to positive values 
of £2 near —E2, and this simplifies the maps we derive below. An analogous 
procedure is used to define local coordinates near ±£3. 

To construct local coordinates near P, we select a cross-section transverse 
to P, say 9\ = 0. Near P, the flow induces a map from that section to itself, 
with P corresponding to a fixed point of the map. We define £1 = r\ — f±, 
where f\ is the value of r\ at the fixed point; £1 is the analogue in the map to 
the radial coordinate for the flow near P. The remaining local coordinates on 
the cross-section are defined by restricting the expanding and contracting di- 
rections at P, as defined above, to the cross-section. Local coordinates can be 
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extended to a neighbourhood of the whole of P in the fully symmetric case by 
applying equivariance under k±. Finally, small symmetry-breaking perturba- 
tions will not change the local structure near P, and we can extend to slightly 
altered local coordinates (£1, #i, X2, £3) in a neighbourhood of P so long as we 
remember that symmetry-breaking terms may have a different effect at each 
value of 0i, so for instance, z\ = fie 101 where f\ = r\{0{) is a function of 0\. 
Note that the global polar coordinates (n , 0\ ) are well-defined near P even in 
the presence of small symmetry breaking since P is far from the origin. 
Cross-sections in M 4 are defined in terms of local coordinates as follows: 



Hf 1 


= {{tl,0l,X2,X 3 


): 


161 < h,0 < 0i < 2tt, \x 2 \ 


< h, \x 3 \ = h} , 




= {(£i,8i,x 2 ,x 3 ) : 


161 < h,0 < 0i < 2tt, \x 2 \ 


= h, \x 3 \ < h} , 




= {(zi,&,x 3 ) : 


z\ 


= h, |6| < h, \x 3 \ < h} , 






= {(21,6,2:3) : 


z\ 


< h, |£ 2 | < h, \x 3 \ = h} , 






= {(21,2:2, 6) : 


Z\ 


< h, \x 2 \ = h, |£ 3 | < h} , 




trout 
^3 


= {(zi, x 2 , £3) : 


Z\ 


= h, \x 2 \ < h, |£ 3 | < h} . 





The cross-sections H 2 n and H™ 1 (resp. and H^ ut ) work equally well near 
±£2 (resp. ±E 3 ) so long as the local coordinate £2 (resp. £ 3 ) is interpreted 
correctly, as described above. 
We also define a Poincare section for the periodic orbit P: 

H? = {(£1, 01, 2:2,2:3) : 161 <h,e 1 =0,\x 2 \ <h,\x 3 \ <h}. 

Trajectories visiting P first cross Hf 1 , may then cross several times, and 
eventually leave the neighbourhood of P on crossing H° nt . 

7.2 Local maps 

Within a neighbourhood of each of ±E 2 , i-^3 an d P, so long as certain non- 
resonance conditions on the eigenvalues are satisfied, the dynamics can be 
linearised using the Hartman-Grobman theorem [37]. In the fully symmetric 
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case, the dynamics near P can be approximated by: 



37 



£1 = -2£, 



1, 



X2 = e±x 2 , 



X3 



where e\ and c\ are positive constants. Without loss of generality, we have 
assumed that the radial eigenvalue is —2, and that the angular speed is 1. 
Solving these equations, we find the local map <pi : Hf 1 — > H° ut is given by: 



X2 

h 



7i 



h In 

ei 



X2 
h 



, /isgn(x 2 ),/isgn(x 3 ) 



X2 

h 



Si 



(13) 

where the initial value of x 3 satisfies \x 3 \ = h, where sgn(x) = +1 if x > 0, 
sgn(x) = —1 if x < 0, and sgn(0) = 0, and where 8± = ci/ei, 71 = 2/ei. 

The argument that symmetry-breaking does not affect this local map goes 
as follows. The transition from Hf 1 — > //° ut has three parts. First, the trajec- 
tory travels from Hf 1 — > i?/ 3 in less than one circuit around P. The trajectory 
does not get very close to P in this time, having started at least a distance h 
from it. Since the Cj's, which control the symmetry breaking, are assumed to 
be much smaller than h, the fully symmetric flow yields an adequate approx- 
imation of the true flow. Second, the trajectory makes ri\ circuits around the 
periodic orbit from Hf to Hf, where n.\ is a non-negative integer no greater 
than T\/2tt. These circuits are governed by the linearised Poincare map and its 
Floquet multipliers: e _47r , e 27rei and e~ 2nCl in the radial, expanding and con- 
tracting directions, respectively, where, to leading order in the ej's, the period 
of P is 2tt. The number m is unchanged by the weakly broken symmetry, and 
so, to leading order, this part of the map is unchanged. Third, the trajectory 
travels from Hf — ► H° nt in less than one circuit around P and again is not 
too close to P, so the fully symmetric flow yields an adequate approximate of 
the true flow. Composing these three parts yields (13), to leading order. 



Local maps ■ H™ 



H% nt and </> 3 : H 3 n -> H$ ut are obtained similarly: 
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where Cj and ej are the absolute values of the real part of the contracting and 
expanding eigenvalues at +Ei, <5j = Ci/ei, ji = 2/ej, and \x 2 \ = h. As for <fii, 
the radial eigenvalues and the angular speeds are chosen to be —2 and 1. 

7.3 Global maps 

The global map ^\ 2 : H° nt — > H 2 n takes orbits from a neighbourhood of P to 
a neighbourhood of +E 2 . We write 

^12(6, 0i,x 2 = h,x 3 ) = (fi = h,9 1 ,£ 2 ,X3) 

and initially do not include symmetry-breaking effects. The unstable manifold 
of P is two-dimensional and, locally, intersects H° nt at 

W u (P)nH° ut = {^ 1 ,9 1 ,x 2 ,x 3 ) : £1 = 0,0 < 6 1 < 2ir,x 2 = ±h,x 3 = 0} (16) 

The manifold W U (P) has two branches: the positive branch intersects H° ut 
with x 2 = h and the negative branch intersects H° nt with x 2 = —h. The 
positive branch forms a connection from P to +E 2 and is the solution we now 
linearise about, while the negative branch forms a connection from P to — E 2 
and will be discussed later. The positive branch of W U (P) intersects at 

{(ri,0i,&,x 3 ) :r 1 = h,0<9 1 < 2tt,&=&,x 3 = 0} (17) 

where £2 is a small constant. The K\ symmetry forces the heteroclinic orbit 
corresponding to the choice 9\ in (16) to have an angular component in H 2 n 
of 6*i + 9\ for some constant 9±, i.e., the global map acts on the angle as a 
rigid rotation. Furthermore, trajectories that are near but not on the unstable 
manifold of W U (P) have £2 and x 3 depending on the initial £1 and X3 but 
not on 9\, while 9\ = 9\ + 9\ where 9\ is a function of the initial £1 and x 3 . 
Equivariance under k 3 ensures that the subspace x 3 = is invariant, that x 3 
is an odd function of x 3 , and that 9\ and £2 are even functions of x 3 . (The 
k 2 symmetry has no role in determining the form of ^12 although it can be 
used to construct a map from P to — E 2 once ^12 is known.) Writing a Taylor 



May 7, 2008 



13:39 Dynamical Systems An International Journal kr 



Effect of symmetry breaking 



39 



series in the small quantities £1 and X3 therefore yields 



6>i + 0i (0,0) +h.o.t. 



&(£l,X 3 ) 



6(o,o) + h.o.t, 

|^(0,0)x 3 + h.o.t. 



x 3 (^i,x 3 ) 



dx 3 



where h.o.t. denotes higher order terms. Effectively, so long as 6\ and £2 are 
non-zero, they can be replaced by constants, while X3 depends linearly on X3. 
We write A\ = §fK0,0) and Bi = 6(0,0), and note that A 1 > since the 
region of phase space with X3 > is dynamically invariant. 

The effect of weak symmetry breaking on these expressions is as follows. 
First, the symmetry x 3 — > —x 3 is broken by including terms that are odd 
in x 3 in the expressions for 6\ and £2, and terms that are even in x 3 in the 
expression for x 3 . We multiply all such terms by an overall factor of £3, which is 
a real constant that controls the magnitude of the breaking of the k 3 symmetry. 
Then the lowest order contribution to 9\ and £2 will be a term in £3X3 while 
x 3 will pick up a term linear in £3. At leading order all quadratic terms can 
be dropped, so the only new term is one linear in £3 in the expression for x 3 . 
Second, breaking the k\ symmetry will result in a weak dependence of all the 
coefficients on 0i, with the dependence being periodic in that variable. We 
introduce the parameter £1, which is a real constant that multiplies all terms 
that break the K\ symmetry and that controls the magnitude of the symmetry- 
breaking terms. For example, A\ will become A\ + £i/i(0), with the caveat 
that this term must remain positive, for all and £1. Third, weakly breaking 
the symmetry X2 — ► —X2 will not affect the form of this map. 

Putting all this together results in a map ^12 : H° nt — > H^ 1 '- 

*i2(fi,0i, x 2 = h,x 3 ) = (n = h,0! = 0i + $1, £2 = B 1: 

X 3 = AiX 3 + £3 + £1X3/1 (0i) + £l£ 3 5'l(6 l l)), (18) 

where $1, A±, B± are constants, and f±, g\ are 2-/r-periodic functions of 0i. 
The 0i dependence cannot be treated using Taylor series expansions, as 0i is 
not a small quantity. We explain below why some quadratic terms (£1X3 and 
£i£3) need to be kept. 
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Similarly, a map from P to — E<i can be constructed. This has precisely the 
form of (18), except that it starts from X2 = —h. Breaking of the k<i sym- 
metry means coefficients in the map will be slightly different but the map is 
unchanged at leading order. 

The map ^23 : — > is calculated in a similar way. In the fully 

symmetric case, we linearise about W n (+E2), which intersects H% nt at (z\ = 
0,^2 = 0, X3 = h) and H™ at (z\ = 0,^2 = ^£3 = £3) where £3 is a small 
constant. For orbits near W u (+E2), the value of £2 at does not influence 
the final position to leading order and z\ at H^ 1 depends linearly on the values 
of z\ at i?2 Ut : z\ = A2^ 2 z\ for real constants A2 > 0, <&2- If the K\ symmetry 
is broken, W u (+E2) leaves H™ 1 with z\ = and arrives at with z\ = ei, 
where <L\ = ei(a r + ia;) for a r and a\ real constants and e\ as defined earlier. 
Writing the resulting map in terms of the real and imaginary parts of z\\ 

*23(n, ^1,6)^3 = h) = (xi = £ia T + A 2 r 1 cos(9 1 + $ 2 ), 

yi = eiai + A 2 n sin(6»i + $ 2 ), 
X2 = /i,| 3 = B 2 ), (19) 

where a r , a;, A2, B2 and $2 are real constants determined by the global flow, 
and A2 > 0. As in ^12, there are 27r-periodic functions of 6\ in the map, 
but here the functions are known explicitly because the z\ variable is small 
throughout the transition from +E2 to +-E3, and the dynamics of z\ is well- 
approximated by a scaled rotation. Similar maps can be obtained for the three 
connections —E2 — ► +E% and ±£^2 — ¥ —E%\ although the coefficients will be 
slightly different in each case, to lowest order we obtain the same map for each 
of the other connections so long as the signs of the X3 (resp. x 2 ) components 
are chosen appropriately on the incoming (resp. outgoing) cross-sections (for 
example, the map from — E2 to — £3 will have X3 = —h and X2 = —h). 
The global map Vl/31 : H% nt — > Hf 1 is calculated in a similar way: 



*3i(ri = Mi,x 2 ,6) = (ll = -B3,£l = 01 + $3, 

x 2 = A 3 x 2 + e 2 + €1X2/3(^1) + eie 2 cfe(0i), 

x 3 = h), (20) 
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where A3, B3 and $3 are real constants, f$ and g% are 27r-periodic functions 
of #1, and €\ controls the size of the terms that break the K\ symmetry. The 
parameter €2 introduced in (20) is analogous to £3, and is a real quantity that 
controls the size of all terms that break the K2 symmetry. Similarly to the case 
for A\ argued above, we take ^3 + ei/3(6'i) to be positive for all values of e±. 
A similar map can be defined near the connections from —E3 to P, and will, 
to leading order, be identical to (20) so long as X3 = h is replaced by X3 = — h. 

The effect of each of the global maps defined above is, at leading order, to 
rotate the angular variable by an order one amount that is independent of 
other variables, to set the radial variable to a constant, and, in the absence 
of symmetry-breaking, to scale the variable that measures proximity to the 
cycle. Symmetry-breaking enters in two ways. First, it destroys the invariant 
subspaces thus destroying some of the heteroclinic connections that made up 
the cycle. Second, breaking the n\ rotation symmetry allows 9\ dependence to 
enter into the maps, most importantly through the variables X3 in the ^12 map 
and X2 in the ^31 map. It is this 9\ dependence that allows the heteroclinic 
tangencies discussed in Section 2. 

7.4 Return maps 

Return maps approximating the dynamics near the heteroclinic cycle can now 
be computed by composing the local and global maps in an appropriate order. 
For instance, to obtain the various forms of the map R : H" 1 — ► given by 
equations (1-5) we calculate R = V>23 4>2 1P12 °4>i° Y^i (p3 i n the usual way. 
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